The Quantum-Classical Correspondence in Polygonal Billiards 
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We show that wave functions in planar rational polygonal billiards (all angles rationally related 
to tt) can be expanded in a basis of quasi-stationary and spatially regular states. Unlike the energy 
eigenstates, these states are directly related to the classical invariant surfaces in the semiclassical 
limit. This is illustrated for the barrier billiard. We expect that these states are also present in 
integrable billiards with point scatterers or magnetic flux lines. 



PACS numbers: 03.65.Sq, 03.65.-w, 05.45. Mt 
I. INTRODUCTION 

The relation between wave motion in the short-wave- 
length regime and the corresponding ray dynamics is of 
fundamental importance in quantum mechanics, electro- 
magnetics and acoustics. It has been found in quantum 
mechanics that the properties of the stationary solutions 
of the Schrodinger equation, the eigenstates \Ej) of the 
Hamiltonian, reflect the degree of order in the classi- 
cal ray dynamics. For example, energy eigenfunctions 
{x,y\Ej) of classically integrable systems with quasiperi- 
odic motion on invariant tori are regular, while eigenfunc- 
tions of classically chaotic systems with ergodic motion 
on energy surfaces are typically irregular. For generic 
systems, it has been conjectured |l| (see ||^ for a review) 
that in the (semi-)classical limit ?i — > the averaged 
Wigner transforms of typical energy eigenfunctions "con- 
dense" uniformly onto the underlying classical stationary 
objects in phase space, which are usually invariant tori, 
chaotic components or entire energy surfaces. In this ar- 
ticle, however, we consider a class of systems with exotic 
classical invariant surfaces for which it is not clear a pri- 
ori whether such a condensation scenario exists. 

The classical free motion inside a planar domain Q 
with elastic reflection at the boundary has a constant of 
motion, Hamilton's function H — p^-{-py-\- V{x, y). The 
particle's mass is 1/2, and V{x,y) = if {x,y) E Q and 
oo otherwise. A billiard with polygonal boundary has a 
second constant of motion K(px,Py) if all angles aj be- 
tween sides are rationally related to tt, i.e. aj — mjir/uj, 
where rrij , > are relatively prime integers . How- 
ever, this does not imply integrability in the sense of 
Liouville-Arnol'd The Poisson bracket {H,K} van- 
ishes identically only if mj — 1 for all j = 1, 2, . . ., so only 
rectangles, equilateral triangles, 7r/2, 7r/4, 7r/4-triangles, 
and 7r/2, 7r/3, 7r/6-triangles are integrable. Critical cor- 
ners with rrij > 1 destroy the integrability in a singular 
way: {H, K} is zero everywhere in phase space except at 
a measure-zero set corresponding to the critical corners; 
the phase space is foliated by two-dimensional invariant 
surfaces H,K — constant as in integrable systems, but 
they do not have the topology of tori. The motion is not 
quasiperiodic and characterized as pseudointegrable |^ . 

The quantum mechanical free wave motion with Di- 
richlet boundary conditions on a rational polygon barely 
reflects the classical pseudointegrability if only individual 



eigenfunctions of the Hamiltonian H = H{x,y,p^,Py) 
are considered. A typical eigenfunction is hard to dis- 
tinguish from eigenfunctions in classically chaotic sys- 
tems ||7|-|lO|. It is in general not an eigenfunction of the 
second operator K K{p^,Py) and has nonzero uncer- 
tainty AK = {{K^) — (fC)^)^/^, since the commutator 
[H, K] does not vanish. Only a subtle signature of pseu- 
dointegrability seems to be contained in individual en- 
ergy eigenfunctions, namely in the distribution of zeroes 
of the associated Husimi function The closeness to 
the quantum mechanics of chaotic systems on the one 
side and to the classical dynamics of integrable systems 
on the other side is a first hint of an unusual quantum- 
classical correspondence. 

A further indication pointing in this direction comes 
from an analogy with the metal-insulator transition of 
the Anderson model in three dimensions. Both kinds 
of systems have energy level statistics close to the semi- 
Poisson distribution In the Anderson model, at the 
transition point between extended states with Wigner 
statistics and localized states with Poisson distributed 
energy levels, this intermediate statistics has its origin in 
the multi-fractal character of the wave functions, which 
are neither extended nor localized ||l^. Analogously, a 
typical energy state in a rational polygon is expected to 
be a fractal in momentum space , localized around an 
energy surface but neither extended nor localized within 
this surface; condensation onto a lower-dimensional in- 
variant surface is ruled out. 

As a final indication we will present a numerical study 
on a particular system, the barrier billiard. Instead of 
Wigner transforms in four-dimensional phase space, we 
study distributions in the two-dimensional (H, _ftr)-space, 
where each point represents a classical invariant surface. 
In this context, we define condensation as AH/ (H) 
and AK/ (K) ^ in the semiclassical limit. This ensures 
that mean values of the quantum operators in highly ex- 
cited states can be interpreted as well-defined values of 
the classical constants of motion. Our central issue is 
to demonstrate that even though the energy states pre- 
sumably do not show such a condensation there is an 
alternative basis of states which do so up to classically 
long times. 

The paper is organized as follows. In Sec.|| we con- 
struct the basis for general rational polygons. Its time 
evolution is discussed in Sec. g. In Sec. we illustrate 
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our statements for the barrier biUiard. We conclude with 
a brief summary in Sec. ^ and an Appendix with details 
on the numerical computations. 



II. CONSTRUCTION OF THE BASIS 

In the following, we will introduce states with rela- 
tively small /S.K at the expense of a nonzero but also 
small energy uncertainty A_ff. As for coherent states in 
the harmonic oscillator [Q, we will minimize the prod- 
uct of the uncertainties involved. To do so, we consider 
the uncertainty relation 



^H^K>\\{[H,Kl 



(1) 



It is an easy matter to show that the equality in (|l|) holds 
in a state \L) satisfying 



{H-{H))\L)^ia{K-{K))\L) , 



(2) 



where i"^ ~ —1, a is an arbitrary real number and 
(. . .) = {L\... \L) with (L|L) = 1. \L) has to be a right 
eigenvector of the operator 



H iaK 



(3) 



L does not commute with its adjoint. It therefore does 
not belong to the class of normal operators, which con- 
tains Hermitian and unitary operators as special cases. 
In general, the right eigenvectors of a nonnormal op- 
erator cannot be used as a basis. Instead, right and 
left eigenvectors together form a biorthogonal basis; see, 
e.g., ijl^. In our case, however, {Ei\L\Ej) is a complex- 
symmetric matrix due to the fact that the Hermitian ma- 
trices {Ei\H\Ej) and Kij = {Ei\K\Ej) can be made real. 
Left and right eigenvectors are therefore identical and 
form separately a nonorthogonal basis. Furthermore, it 
can be shown that AH = |a|AAr. In order to have equal 
uncertainties we choose a — 1 (states with a = —I are 
identical). The real and imaginary parts of the complex 
eigenvalues L have a physical interpretation as mean en- 
ergies {H) and —{K), respectively. Note that the clas- 
sical function H — iK is a constant of motion whose 
constant-level surfaces are the invariant surfaces. 

We derive now the properties of the |L)-states in the 
limit h ^ Q. Without loss of generality, we stipulate that 
K and K are homogeneous functions of the momenta of 
degree two, like H and H . Starting from the expansion 
(see, e.g., [HQ) 

^x,Py\i[H,K]\x,y){x,y\p^,Py) = n{H,K} + 0{ri^) , 

(4) 

we will exploit only the fact that the Poisson bracket 
{H,K} is everywhere zero except at isolated critical 
points in position space. The following line of reason- 
ing is therefore also true for integrable billiards with a 
finite number of magnetic flux lines or point scatterers 
(we have checked this for Seba's billiard [18[). In the 



semiclassical regime, both sides of Eq. (^) are very small 
in the region Q\C excluding the union C of neighbor- 
hoods of the critical corners, the area of which shrinks to 
zero as ?i ^ 0. Manipulating Eq. (||) and restricting the 
integration over an |L) -state to region Q\C gives 



(5) 



with {AHQ\cf - {{H - {H)) 



Q\Ci 



etc (Note that 



1)q\c < !)■ The smallness of the l.h.s. of Eq. (^) carries 
over to both sides of Eq. (||). From this we conclude that 
a function {x,y\L) in region Q \ C is locally either very 
small or can be approximated by a joint eigenfunction 
of both operators H and K. Such a joint eigenfunction 
cannot fulfill the boundary conditions globally, otherwise 
the billiard would be integrable. Loosely speaking, re- 
gion Q\C must be divided into subregions in each of 
which a joint eigenfunction (or a vanishing function) ful- 
fills the boundary conditions locally. Joint eigenfunctions 
of neighboring subregions match smoothly, so they must 
have roughly the same number of nodal lines. Hence, 
|L)-states have features of energy states of integrable 
systems: (i) the functions {x,y\L) are in some regions 
regular while in other regions, separated by "caustics", 
vanishing; (ii) they can be labelled by two "quantum 
numbers" (^1,722); (iii) the eigenvalues L = (H) — i{K) 
are regularly distributed in the complex plane. The last 
property holds because {H) (all arguments are also valid 
for quantities derived from K) is asymptotically equal to 
(H) Q\c I which is approximately a homogeneous function 
of the quantum numbers of degree two due to the homo- 
geneity of H. The non-relevance of {H)c can be un- 
derstood with a renormalization procedure. Going from 
(77.1,77.2) to (2ni,2?i2), we get a new wave function with 
slightly larger region Q\C containing four times more 
nodal lines, and region C being a four times smaller copy 
of the old region C. Hence, (iif)g\c roughly increases by 
a factor of four, while {H)c stays constant. This leads to 
(H) —^ (H) as more and more renormalization steps 
bring us towards the semiclassical limit. 

The approximate homogeneity of (H) q\c with respect 
to the quantum numbers ensures that the local transi- 
tions between different joint eigenfunctions induce only 
a small uncertainty {AHq\c)^ of order (H) ^ {H'^), 
increasing roughly by a factor of four under the ac- 
tion of the renormalization. The same factor is an up- 
per bound for the increase of (AHc)^, corresponding to 
a (Jf^)-scaling weighted with the size of C. The sum 
{AHY — {AHq\cY + {^HcY therefore increases by a 
factor of four. Hence, 



{AHf oc {H) , {AKf cx {K) 



in the semiclassical limit. Consequently, 



AH „ 
lim T^r^T * , lim — > , 



h^o (H) 



h^o IK) 



(6) 



(7) 



i.e. the |L)-states condense onto the invariant surfaces. 
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III. TIME EVOLUTION 

The time dependence of |i)-states is non-trivial since 
L does not commute with the Hamiltonian; \L{t)) is in 
general not an eigenstate of i for i > 0. Let us define 
three time scales associated with a given state \L) at time 
t = Q: the quantum mean period, the lifetime and the 
classical mean free time 



2-1:% 

Th) 



2'Kh 
'ah 




(8) 



where A is the area of the billiard. From Eq. (|6|) we get 
tq ^ Tl ~ Tc, i.e. the state is quasi-stationary with a 
lifetime of order of the classical mean free time. 

At first glance, it seems that the state fails to con- 
dense onto an invariant surface for classical long times 
t ^ tq. However, the requirements for condensation, 
small relative uncertainties and constant mean values of 
H and K, may persist beyond the lifetime of the state. 
Clearly, this is true for the relative uncertainty and the 
mean value of H since this operator commutes with the 
evolution operator exp (—iHt/h). Keeping in mind that 
{i[H, K]) is relatively small, we may expect that the rel- 
ative uncertainty and the mean value of K change slower 
than the state itself. To see this, we define the time scale 
associated with 

{K){t) = {L{t)\K\L{t)) 

OC 

= ^ {L\E,)K,,e^^^^-^''^'^\E,\L) . (9) 

in the same way as the lifetime in Eq. as tk — 
2TTfi/AHii, where 

is the mean energy difference, i.e. mean frequency dif- 
ference, of the energy states involved, but in contrast 
to AH the states are weighted according to their rel- 
evance for {K){t). Two extreme cases show that defi- 
nition ( p^ ) is reasonable: a diagonal matrix Kjk gives 
AHk = whereas a uniform matrix Kjk ~ constant 
gives AHk = AH. 

The series (^) can be expanded in orders of t, with t 
much smaller than tk but with the possibility of being 
larger than the quantum mean period and even the life- 
time. Clearly, the zeroth-order term is {K). The next 
terms are roughly of order {K)AHKt, {K){AHKt)'^ , and 
so on. If we compare this to 



{K){t) = {K) + m,K])t/h + 0{t') 



(11) 



we find that the order of tk can be estimated as 
2tt (K) / {i[H , K]) . It follows that in the semiclassical 
limit (K) {t) does not change relative to its initial value 
for classically long times t below tk S> tc, even though 
\L{t)) may differ strongly from the initial state \L). 



Following the same line of reasoning shows that the 
time below which {K'^)(t) relative to its initial value 
is constant scales as {K'^) / {i[H , K'^]) . Starting from 
Eq. (^) the following equation can be derived 

m.K^])-m)m.K])^2{{K-{K)f) . (12) 

With the renormalization procedure it can be shown that 
the asymptotic behavior of the r.h.s of this equation is 
bounded from above by {K) ^ . Together with the already 
known result {i[H,K]) cx (K) we finally get the upper 
bound {i[H,K'^]) cx {K'^). Hence, the time scale which 
governs {K'^){t) is of the same or of larger order as tk. 
From this we can conclude that for times smaller than tk, 
AK(t) remains small if compared with {K){t). Hence, 
Eq. (||) holds not only for times of order tl ~ tq but also 
for times much larger than tq and well below tk- The 
condensation of |L)-states onto invariant surfaces outlives 
the lifetime of the states up to classically long times. 



IV. EXAMPLE: THE BARRIER BILLIARD 

We illustrate all statements for the barrier bil- 
liard [p9|-pl|, a rectangle with width and height ly 
and a vertical barrier of length ly/2 placed on the sym- 
metry line X = lx/2] see Fig. 0(a). The billiard is not 
only pseudointegrable, it is also almost-integrable [^ , 
i.e. is composed of several copies of a single integrable 
sub-billiard, here the rectangle shown in Fig. 0(b). The 
function K = is a second constant of motion. The 
general formula for the genus of the invariant surfaces 
gives 2, i.e. the surfaces have the topology of two-handled 
spheres and not that of tori (single-handled spheres). 



(a) 



(b) 



Ix 



region 2 



region 1 



L/2- 



FIG. 1. (a) Barrier billiard, rectangle with a barrier between the 
points {x,y) = {lx/2,0) and (Ix /2, ly /2). (b) Symmetry reduced 
system. The nontrivial wave functions fulfill Dirichlet (Neumann) 
boundary conditions on solid (dashed) lines. 

The energy eigenfunctions are solutions of the Helm- 
holtz equation with Dirichlet boundary conditions on 
the polygon. The functions are odd or even with re- 
spect to the symmetry line. The odd ones are trivial 
eigenfunctions of the integrable sub-billiard. We there- 
fore deal only with the even ones, which fulfill mixed 
boundary conditions on the symmetry reduced polygon; 
see Fig. 0(b). We have calculated the solutions numeri- 
cally with the mode-matching method for the parameters 



1, Ix = ttvott/S, and ly = SvStt/tt as described in 



App. Al. The statistical properties of the energy levels 
are found to be close to the semi-Poisson distribution [n2| . 
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The energy eigenfunctions are not eigenfunctions of the 
operator K = p^. The numerical computation of the un- 



certainty AK is explained in App. A 2, Figure g shows 
that there is no trend towards vanishing relative uncer- 
tainties. This indicates that energy states of the barrier 
billiard do not condense onto the invariant surfaces. 
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FIG. 2. Relative uncertainty AK/{K) of energy states \Ej). 
A local Gaussian average with a variance of 200 is performed. In- 
set: Contour plot of the probability density associated to the 626th 
energy state with E = 646.03, (K) = 135.05 and AK = 144.01. 



apart from a localization around the critical corner, the 
type-B functions are restricted to the lower (ni odd) or 
upper {ni even) half of the billiard, bounded by a caustic- 
like curve. 
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FIG. 3. Eigenvalues L = (H) — transformed into the ac- 

tion space of the sub- billiard according to Eq. (p^. The dotted line 
marks the centre of the transition region. The solid lines indicate 
parts of the EBK lattices defined in Eqs. (|l^-(pi|). 



The nontrivial |L) -states are also calculated with the 
mode-matching method. As explained in App. A3, we 
cannot compute as many of these states as energy states. 
However, even in the accessible low-energy regime our 
theoretical results on the semiclassical behavior of these 
states will be well confirmed in the following. 

The regular pattern of the eigenvalues L = (H) — i{K) 
can be most clearly seen when transformed into the real 
action variables of the integrable sub-billiard 



I. 



ZTT TT 



(K] 



(13) 



Note that the topology of the invariant surfaces in the 
full system rules out action-angle variables As can 
be seen from Fig. || the "action space" is split into two 
regions A and B separated by a transition region. Away 
from this transition region, and the lines — and 
ly = 0, the eigenvalues are approximately located on 
regular lattices which are given by EBK-like quantiza- 
tion rules: 



ilx,ly) = ni 



1 n 



in region A and 



ni 1 „ 7\ , 

— + -,2^2 + - U 

2 2' 4/ 



(14) 



(15) 



in region B with ni,n2=0,l,.... For example, we have 
{Ix,Iy) ~ (6.75,37)?i for the eigenfunction of type A 
and {Ix,Iy) ~ (21,9.7)?i for the function of type B in 
Fig. ^. Both types have relatively small uncertainties 
(see App. A 4 for numerical details) and look rather reg- 
ular. While type- A functions cover the billiard uniformly. 



^ as) !^ <ss> 1^ <^ 3^ <^ <^ g 

SS) ^ ^ ^ 

_ - =. =: 

— ^ ^ 3^ ^ ~ ^ -5 ^ ~" ~ 
iZ. rzi. = ^^^^^^-cL^ — 

■^^ ^ "5^ 5" ,S S ^ S ^ '3> 


iiiiiiiiiiili^BMiiiiHwiiiPHii 






S — ^ i 3. $ 

S, S J B S g ^ 


S J 8 g g S S 

^ ^ <sS) <^ !0> 


iiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiM^^^ 
iiiiiiiiiiiiiiiiiiite 
iiiiiiiiiiiiiiiiiM^^^^^ 
iiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiM^^^^^ 







FIG. 4. Probability density of a state |L) of type A (left) with 
eigenvalue 656.30 - i64.93 and AH = AK = 34.88; type B (right) 
with eigenvalue 670.95 - i630.52 and AH = AK = 25.14. 

The eigenvalue pattern in Fig. || resembles strongly 
that of integrable systems with separatrices [p3|-[25|. It 
is therefore not surprising that the statistical properties 
of the mean energies of |L)-states are similar to those of 
energy levels of integrable systems. For example. Fig. |^ 
confirms that the nearest-neighbor statistics is in agree- 
ment with the Poisson distribution. 
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FIG. 5. Probability density P{s) of the spacing s between ad- 
jacent values of the first 800 mean energies (H). The data is well 
fitted by the Poisson distribution (dotted line). 

The scaling laws ^ for the uncertainties are verified 
in the following way. First note that lines of constant 
mean energy (ff) are circles in the scaled action space 
(2Ix/lx, ly/ly)- This plane can therefore be conveniently 
parametrized by the radial coordinate (H) and the polar 
angle cf) = arctan [(/j,/^j,)/(2/3;//2;)]. The first scaling law 
holds if (AH)'^ divided by (H) is a function of alone. 
Figure |^ shows that this condition is well satisfied. This 
confirms also the second scaling law since AK = AH 
and (K) scales as (H). 
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FIG. 8. State of type B with L = 670.95 - i630.52, cf. Fig 
at time t = lOrq (left) and i = tl 26.7tq (right). 
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FIG. 6. (/\Hf/{H) vs. 4> with 400 < (H) < 800. The maxi- 
mum lies in the transition region of action space; cf. Fig. ^ 



\L) 



The main features of the time dependence of 
states can be observed in Figs. and ^; see App. A_5 
for the numerical aspects. For small times t around 
TQ <C TL « 26.7rQ, the time evolved state \L{t)) and 
the initial state \L) are similar; the overlap |(L(<)|L)p is 
close to one. For times < ~ tl both states differ consider- 
ably; the overlap is close to zero. The state has lost some 
of its regularity, but {K){t) and AK{t) stay constant up 
to order {K) for longer times well below tk ~ 2736tq. 
Interestingly, the states of type B become more uniformly 
distributed in configuration space in the course of time, 
cf. Figs. |(b),|(a) and (b). 




FIG. 7. Overlap \{L(t)\L)\^ (solid), {K){t)/(K) (dotted) 
and AK{t)/{K){t) (dashed) for the state of type B with 
L = 670.95 - j630.52; cf. Fig. Q. 



V. CONCLUSION 



We have formulated the quantum-classical correspon- 
dence in rational polygonal billiards in the following way: 
there exists a basis of quantum states with each state 
condensing individually onto a classical invariant surface 
H,K — constant up to classically long times in the sense 
that the relative uncertainties AH/{H) and AK/{K) 
vanish in the semiclassical limit. We have presented some 
hints, like the analogy to the Anderson metal-insulator 
transition, and numerical evidence for a particular sys- 
tem, the barrier billiard, that the energy states do not 
show such a condensation. We have then introduced an 
alternative basis of states for which we have explicitly 
shown that (i) they are quasi-stationary; (ii) they con- 
dense onto the invariant surfaces up to classically long 
times even exceeding their lifetimes; (iii) in configuration 
space they show a regular nodal-line structure, possibly 
with caustic-like boundaries; (iv) the eigenvalues form a 
regular pattern in the complex plane. 

Whether these states condense uniformly onto the in- 
variant surfaces in phase space is another question. At 
first sight, our numerical results on the barrier billiard 
seem to indicate that there cannot be uniform condensa- 
tion since a fraction of the states cover only one half of the 
configuration space. However, each such a state becomes 
more uniformly distributed on the invariant surface for 
times beyond the lifetime but remains in the neighbor- 
hood of the surface for classically long times, i.e. there 
can be uniform condensation after some transient time. 

It is important to note that the condensation scenario 
holds for general rational polygonal billiards and also if 
the time-dependent Schrodinger equation is replaced by 
the wave equation used in acoustics and electromagnet- 
ics, even though the individual time scales are different. 

I would like to thank M. Sieber and J. Nockel for dis- 
cussions. 



APPENDIX: NUMERICAL COMPUTATIONS ON 
THE BARRIER BILLIARD 
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1. Energy states 

We compute the energy eigenvalues and eigenfunctions 
with the mode-matching method; see, e.g., ||^, p6| , p7t . Let 
us first set h = 1 and then divide the symmetry-reduced 
barrier bilUard in two regions as shown in Fig. ^b): re- 
gion 1 with y < I and region 2 with y > I. The length of 
the barrier, I, is in our case fixed to ly/2 but the following 
derivations hold also for general < I < ly. In region 1 
we expand the wave function as 

oo 

$1 = ^ am sin {2rmrx/lj;) sin {g2my) (Al) 



m— 1 



with 



(A2) 



This kind of expansion is nontrivial since gj can be 
imaginary and a„j complex. By construction, the func- 



tion (Al) fulfills the Helmholtz equation -V^$i = E^i 



with Dirichlet boundary condition on x = 0, y = 0, and 
X — lx/2. In region 2 we take an analog function, 

oo 

$2 = ^ bm sin [(2m - l)Trx/lj;] sin [g2m-i{ly - y)] , 

m— 1 

(A3) 

which satisfies Dirichlet boundary condition on x = and 
y = ly but Neumann boundary condition on a; = lx/2. 
We stipulate that both functions match smoothly at 
y = I, i.e. for < x < lx/2 we require 



$i(a;,0 = $2(a;,0 



and 



dy 



(x,l) 



d^2 
dy 



{x,l) 



Inserting the identity 

oo 

sin {2nz) — Ank sin [{2k — l)z\ 



(A4) 



(A5) 



(A6) 



k=l 



for < z < 7r/2 with the orthogonal matrix 

2 / sin [{2n - 2k + l)7r/2] sin [{2n + 2k - l)7r/2] 



A. 



2n - 2fc + 1 



2n + 2fc - 1 

(A7) 



into the first matching condition (A4) and solving for the 
coefficients of sin [{2n — l)Trx/lx], n = 1, . . ., gives 

oo 

^ am sin {g2ml)^mn = K sin [g2n-l{ly - 0] • (A8) 

m— 1 

Similarly, we get from the second matching condi- 
tion (|A5|) 



a™g2m cos {g2ral)Aran = -^iff2n-l COS [g2n-l{ly " 0]' 

rn—1 

(A9) 



We now rewrite the relations (A8)-( A9) by using the def- 
initions 



am = amg2m COS {g2ml) , 

bm = bmg2m-l COS [g2m-l(ly - 0] 



and the real function 



as 



tan (g„J) 

5m 



(AlO) 
(All) 



(A12) 



0^mf2m{E, l)Amn = inf2n-l{E, ly — I) (A13) 



and 

oo 
m— 1 

The last two equations are combined to 

oo 
m— 1 

with the real matrix 

MmniE) = [/2™(^,0 + /2n-l(S,;y-0]A 



(A14) 



(A15) 



(A16) 



Equation (A_15) has a solution provided det M(i?) = 0. 
Finding the energy eigenvalues Ei is therefore equivalent 
to finding the zeroes of the determinant of M . We ap- 
proximate M by a 500 x 500 matrix which is sufficient 
for calculating the first 100 000 zeroes. Figure ^ shows 
the determinant as a function of the energy. It is conve- 
nient to search roots numerically only between two con- 
secutive poles. We therefore rewrite Eq. ( A12 ) using the 
identity 



- — tan — - ^ 



4z 



^ Z2 _ (On - 1)2 
n— 1 ^ ^ 



as 



2 °° 1 

n=l ^ I. ' ^21 ) 



(A17) 



(A18) 



where the summation is only over odd n. From this ex- 
pression the poles can be easily read off. Note that the 
poles are not degenerate if l^/ll and {ly — O^/'k ^'^'c irra- 
tional numbers. The interval between two given consecu- 
tive poles is divided into 400 subintervals, each assumed 
to contain at most a single change of sign of det M. The 
bisection method is then employed in order to find each 
zero with an accuracy 10~^ of the mean level spacing. 
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FIG. 9. Determinant of matrix M in Eq. (A.l£) 



Having determined the energy levels, we get for each of 
them the real quantities hm and hm from Eqs. (A14) and 
( A15| ) and then the wavefunction ( |A1[ ) and (A3), which 
we finally normalize to unity for the full billiard. 



2. Uncertainties of energy states 

We here determine the uncertainty (Aii')^ = {K'^) — 
(K)^ in a given \Ej)-state. We take $i = ^jjEj), 
g = g{Ej), and $2 = '^2iEj) from Eqs. (^)-(^ to 
compute Kjj — {Ej\K\Ej) = (K) from 



1^/2 fly / 02 



/o Ji \ dx 
A straightforward calculation gives 



Tr^$2 dxdy 



m— 1 

00 



sin (252mO 



2ff2 



(A19) 



(A20) 



sin [2ff2m-l(^y - Q] 
252m-l 



compute L in energy-state representation, i.e. we calcu- 
late the matrix elements {Ej\L\Ek) = EjSjk —iKjk with 
Kjk = {Ej\K\Ek) given by Eq. ( |^20| ) if j = k, otherwise 



K 



jk 



/ ^ UmUmil'2m 



sin (g2mO sin (.g+„0 



52r, 



52r, 



J X! blnbmh2m~l X 



(A22) 



sin [ff2m-i('!/ - 0] sin [5^-1(^1/ " 0] 



52m-l 



52m-l 



with at„ = a,„ = am{Ek), = KniEj), Ki = 

bmiEk), .g+ = g*^{Ej) + g^{Ek), 5m = 9*^{Ej)- grn{Ek), 
and /im = (mTr/Zj.)^. A LAPACK routine is used to diag- 
onalize the complex matrix (E'j with j^k < 1000, 
giving a rough approximation to the first 800 eigenvalues 
of L. We use them as initial guesses for Newton's method 
in order to find reliable approximations to the eigenval- 
ues. Finally, we compute for each determined eigen value 
t he c omplex quantities fim and 6,„ from Eqs. ( A14|) and 
( A15| ) and then the normalized wavefunction (Al) and 



4. Uncertainties of |L)-states 

The main advantage of using the mode-matching 
method for the |L)-states is that we get the uncertainties 
AK and AH as accurately as for the energy states in 
App. We take <I>i = <l>i(L,-), $2 = ^{Lj), and 

g = g{Lj) from Eqs. @, (|a|), and ( |A211 ) to com- 
pute firstly (K) starting from Eq. ( A19| ). A straight- 
forward calculation shows that (K) is given by the 
r.h.s. of Eq. ( |A22D with aj^ = a*„{L.j), Um = ara{Lj), 
= b*^iEj), Ki = 6m (ij), gti = 9*n{Ej) + ffm(ij), 
5m = 9m{Ej) - 9m{Lj), and hm = {rmr/l^)'^. Simi- 
larly, we get for {K'^) = (p^) the same equation but with 

and for (H'^) we use h„ 



{rmr/lx)^; for {H) we use 

\9fniEj 



■ 9m{Lj) + {rmr/l^y 



with am = a,n{Ej), hm = hm{Ej), gm = gm{Ej), and 
hm = {rmi/lxY- Analogously, we get for {K"^) = (p!^) 
again Eq. ( |A20| ) but with hm = {rm: /Ix)^- 



3. |L)-states 

Eigenvalues and eigenfunctions of L can be co mpu ted 
with the procedure described in App. Al if Eq. (A2) is 
replaced by 



5^(i)=i-(l-z)(f 



(A21) 



with complex number L. M in Eq. (AI6) is then a 
complex matrix. Finding the complex roots of det M 
is much more cumbersome than finding real roots as in 
the case of the energy states. Because of this, we first 



5. Time dependence of |L)-states 

We here compute {x,y\L{t)), {L{t)\L), {K)(t) and 
{K'^){t). These expressions can be written in terms of 
energy states as 



{x,y\L{t)) 



mm 



00 



{E,\L){^ME, 



Eit 



\{EAL)\' 



(A23) 



(A24) 



and as in Eq. (||). The first 1400 energy states are incor- 
porated in these sums. As in the previous sections it can 
be shown that {Ej\L) is given by the r.h.s. of Eq. ( [A2^ ) 
with al^ = am{E-j), am = am{L), b]^ = Kn{Ej), bm = 
bm{L), g+ = gmiEj) + gm{L), g- = 9m{Ej) ~ 9rn{L), 
and hm = 1. {Ej\K'^\Ek) is also given by the r.h.s. of 
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Eq. ( [A2^ ) if i 7^ fc with at„ = a*^{Ej 



bi 



5m = 9*m.{Ej) - gm{Ek), and K 
one has to use Eq. ( |A2Cl[ ). 



= 9*m{Ej) + gm{Ek), 
{rmr/lx)'^- For j = fc 



[1] 
[2] 

[3] 
[4] 
[5] 



[6] 
[7] 

[8] 
[9] 

[10] 

[11] 
[12] 



M. V. Berry, J. Phys. A 10, 2083 (1977). 

M. Robnik, Nonlinear Phenomena in Complex Systems 

(Minsk) 1, 1 (1998). 

A. Hobson, J. Math. Phys. 16, 2210 (1975). 
A. Zemlyakov and A. Katok, Math. Notes 18, 760 (1976). 
V. I. Arnol'd, Mathematical Methods of Classical Me- 
chanics, Vol. 60 of Graduate Texts in Mathematics 
(Springer, Berlin, 1978). 

P. J. Richens and M. Berry, Physica D2, 495 (1981). 
D. Biswas and S. Jain, Phys. Rev. A 42, 3170 (1990). 
P. Bellomo and T. Uzer, Phys. Rev. A 51, 1669 (1995). 
J. Vega, T. Uzer, and J. Ford, Phys. Rev. E 52, 1490 

(1995) . 

M. Tomiya and N. Yoshinaga, J. Stat. Phys. 83, 215 

(1996) . 

D. Biswas and S. Sinha, Phys. Rev. E 60, 408 (1999). 

E. Bogomolny, U. Gerland, and C. Schmitt, Phys. Rev. 
E 59, R1315 (1999). 



[13] M. Schreiber and H. Grussbach, Phys. Rev. Lett. 67, 607 
(1991). 

[14] R. Glauber, Phys. Rev. 131, 2766 (1963). 

[15] J. Wilkinson, The Algebraic Eigenvalue Problem, Mono- 
graphs on numerical analysis (Clardendon Press, Oxford, 
1965). 

[16] W. Brenig, Statistische Theorie der Wdrme (Springer, 
Berlin, 1992). 

[17] M. Hillery, R. O'Connell, M. Scully, and E. Wigner, 

Physics Reports 106, 121 (1984). 
[18] P. Seba, Phys. Rev. Lett. 64, 1855 (1990). 
[19] R. Zwanzig, J. Stat. Phys. 30, 255 (1983). 
[20] J. H. Hannay and R. J. McCraw, J. Phys. A 23, 887 

(1990). 

[21] J. Wiersig, Phys. Rev. E 62, R21 (2000). 

[22] E. Gutkin, Physica 19D, 311 (1986). 

[23] P. H. Richter, H. R. Dullin, H. Waalkens, and J. Wiersig, 

J. Phys. Chem. 100, 19124 (1996). 
[24] H. Waalkens, J. Wiersig, and H. R. Dullin, Ann. Phys. 

260, 50 (1997). 

[25] H. Waalkens, J. Wiersig, and H. R. Dullin, Ann. Phys. 
276, 64 (1999). 

[26] H. Wu, D. Sprung, and J. Martorell, Phys. Rev. E 51, 
703 (1995). 

[27] P. Exner, P. Seba, M. Tater, and D. Vanek, J. Math. 

Phys. 37, 4867 (1996). 
[28] I. S. Gradshteyn and I. M. Ryzhik, Tables of Integrals, 

Series, and Products (Academic Press, New York, 1965). 



8 



